Out-of-plane equilibrium spin current in a quasi-two-dimensional electron gas under 

in-plane magnetic field 
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Equilibrium spin-current is calculated in a quasi-two-dimensional electron gas with finite thickness 
under in-plane magnetic field and in the presence of Rashba- and Dresselhaus spin-orbit interactions. 
The transverse confinement is modeled by means of a parabolic potential. An orbital effect of 
the in-plane magnetic field is shown to mix a transverse quantized spin-up state with nearest- 
neighboring spin-down states. The out-off-plane component of the equilibrium spin current appears 
to be not zero in the presence of an in-plane magnetic field, provided at least two transverse-quantized 
levels are filled. In the absence of the magnetic field the obtained results coincide with the well- 
known results, yielding cubic dependence of the equilibrium spin current on the spin-orbit coupling 
constants. The persistent spin-current vanishes in the absence of the magnetic field if Rashba- and 
Dresselhaus spin-orbit coefficients, a and P, are equal each other. In-plane magnetic field destroys 
this symmetry, and accumulates a finite spin-current as a — > j3. Magnetic field is shown to change 
strongly the equilibrium current of the in-plane spin components, and gives new contributions to the 
cubic-dependent on spin-orbit constants terms. These new terms depend linearly on the spin-orbit 
constants. 

PACS numbers: 74.78.-w, 74.62.-c, 74.70.Kn, 74.50.+r 



I. INTRODUCTION 



A central goal of spintronics research is an achievement 
of an electron spin manipulation by means of an exter- 
nal electric field [l[ instead of a magnetic field, which 
is widely used now in semiconducting devices for en- 
hancement of an information processing speed. An elec- 
tric field controlled spin-orbital coupling is a promising 
tool [2j in realization of futuristic spin transport de- 
vices. During the last ten years, there has been im- 
pressing progress in both experimental and theoretical 
understanding of the spin dynamics in quantum wells on 
the base of, particularly, narrow-gap semiconductors with 
high g-f actor, and metal-oxide-semiconductor field-effect- 
transistor (MOSFET) structures. Although the dimen- 
sionless spin-orbit (SO) coupling parameter in vacuum 
is as small as Ep / (moc 2 ) ~ 10~ 6 , where Ep ~ leV is 
the Fermi energy of an electron and ttiqc 2 ~ lMeV is 
the Dirac gap, the large value of the SO coupling energy, 
comparable with the Fermi energy, can be ensured by 
large potential gradient on the semiconductor/insulator 
interface of these structures in the presence of macro- 
scopic structural inversion asymmetry (SIA) . Indeed, the 
gate potential applied across the substrate in MOSFET 
results in inhomogeneous space charge distribution near 
the semiconductor/insulator interface. The nonuniform 
macroscopic potential, confining the electrons near the 
interface, varies over a wide range, ~ 10 4- 1000 nm, 
with larger potential gradient, which originates so-called 
Rashba SO interaction [3, 4]. On the other hand, higher 
value of the SO coupling is achieved by choosing special 
semiconducting materials with a bulk inversion asymme- 
try (BIA) in their crystalline structure, where the gradi- 



ent of the crystal potential is large. Most prominent semi- 
conducting compounds have either zinc-blende structure, 
like GaAs and most of III-V compounds, or wurtzite 
structure in II- VI compounds with BIA. Lack of the bulk 
inversion symmetry in these compounds was shown by 
Dresselhaus [5] to originate another macroscopic SO in- 
teraction. 

Effects of both Rashba- and Dresselhaus- SO couplings 
to the physical properties of two-dimensional (2D) elec- 
tron gas are not trivial even in the absence of an exter- 
nal magnetic field. Existence of a SIA or a BIA in a 
disordered 2D system changes reversely the sign of the 
phase-coherent localization correction to the conductiv- 
ity (B-Q, driving the system from a weak localization 
regime into an antilocalization one. Rashba- and Dres- 
selhaus SO interactions equally and independently con- 
tribute to the weak antilocalization correction. Contribu- 
tions of Rashba and Dresselhaus SO interactions to the 
D'yakonov-Perel's spin relaxation rate Q were shown to 
be also additive [loj . On the other hand, the anisotropic 
contribution to the conductivity tensor ll| in the pres- 
ence of both Rashba and Dresselhaus SO terms, the ab- 
sence of spin polarization and suppression of spin accu- 
mulation especially at the equal values of the coupling 
constants a = ±/3, [l2-14|, restoration of the weak lo- 



calization regime back at a = manifest an 

existence of the interference between Rashba and Dres- 
selhaus SO interactions. Although the SO coupling gen- 
erally breaks the spin rotational symmetry, a new type of 
SU(2) symmetry appears [l8| in the case of a = /3, which 
renders the spin lifetime. In the presence of Rashba and 
Dresselhaus terms with equal strength, the SO interac- 
tion rotates electron spins around a single fixed axis. The 
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spin along this axis becomes conserved, nevertheless spin 
aligned in the perpendicular directions undergoes a deter- 
ministic rotation depending only on the initial and final 
points of their trajectory. 

Different experimental techniques have been developed 
recently to control a coupling of spin to the electric field 
[l9l - |2l| . An efficient g-tensor modulation resonance, ob- 
served in a parabolic Al x Ga\^ x As quantum well (l9j with 
varying Al content x = x(z) across the well, provided an 
opportunity to manipulate electron spins by means of 
various electron spin resonance type techniques. An in- 
plane magnetic field in all of these experiments seems to 
be rather favorable for getting a pronounced spin reso- 
nance. SO interactions in a 2D electronic system pro- 
duce an effective in-plane field, which results in an drift- 
driven in-plane spin polarization 22J. An external in- 
plane magnetic field appears to be not always summed 
algebraically with SO induced effective field, and can re- 
sult in the surprising out-of-plane spin polarization 23| . 
which has been observed in a strained n — InGaAs film 
[20j |. On the other hand, Han le p recession of optically 
oriented 2D electrons in GaAs [24j is well described by a 
total in-plane field, given as a sum of the external- and 
SO effective fields. All these facts show nontrivial effects 
of in-plane magnetic field on spin dynamics in quasi-2D 
systems. Effect of an external magnetic field, aligned in 
the normal direction to the electron gas, has been studied 
very well, since the problem can be solved exactly for a 
non-interacting electron gas in the presence of one of the 
SO interaction. In the previous activities, the selective 
coupling of the in-plane magnetic field to the electronic 
spin degree of freedom in the presence of SO interactions 
has been used to probe the interplay of Zeeman split- 
ting with the SO coupling. The in-plane magnetic field 
in non-ideal 2D systems with a finite width, which is 
particularly relevant for heterojunctions and MOSFET 
structures, couples also to the orbital motion, and can 
considerably modify the physics involved. It is there- 
fore important to characterize the various physical effects 
generated in the presence of a parallel magnetic field, in 
order to gain a better understanding of the influence of 
orbital magnetic effects on the physical properties of an 
electron gas with SO interactions. 

Although the electron gas, formed on the semiconduc- 
tor/insulator interface in the heterojunctions and MOS- 
FET structures, has a finite thickness 25j, in the most 
activities concerning the SO interactions it is taken as 
a strictly 2D object by neglecting the finite thickness. 
The thickness of the confined electron gas in these struc- 
tures is varied in the large interval by the gate potential, 
applied across the electron gas, from 10 nm in the inver- 
sion regime up to 1000 nm in the depletion regime. The 
charge distribution and the electron gas thickness can 
be experimentally measured and theoretically estimated 
with high accuracy by means of the self-consistent so- 
lutions of Schrodinger and Poisson equations under the 



charge balance condition. Finite thickness of the electron 
gas was recently suggested by Rashba and Efros (26. 27| 
in order to study the time-dependent gate voltage ma- 
nipulation of electron spins in MOSFETs and quantum 
wells, since the spin response to a perpendicular-to plane 
electric field can be achieved due to a deviation from 
strict 2D limit. 

In this work, we report on our investigation of both 
orbital and spin effects of in-plane magnetic field in a 
quasi-two-dimensional (quasi-2D) electron gas with a fi- 
nite thickness on the spin precession and splitting in the 
presence of Rashba and Dresselhaus SO interactions. In 
this paper we calculate persistent spin current. Note 
that the model has been considered in our previous paper 
H, 2l| in order to study the energy spectrum and the 
Fermi surface under variations of SO coupling constants, 
the gate electric field, the magnetic field and g-factor. 
Generation of a spin flux and its change under exter- 
nal destructive factors is still a controversial issue [3(| 
in spintronics. Generation of a dissipationless transverse 
spin current or a spin Hall current by a driving electric 
field E was predicted 3l|, |32[ in a clean, infinite and ho- 
mogeneous structural inversion asymmetric 2D system. 
Even an arbitrary small concentration of non-magnetic 
impurities was shown [33T4361 ] to suppress totally the uni- 
versal value of the spin Hall conductivity peculiar to a 
clean system. As it was shown by Rashba [33], a spin 
current in the presence of Rashba SO coupling appears 
even at equilibrium in the absence of an external elec- 
tric field, though it does not result in any accumulation 
of spin. A universal equilibrium spin current was shown 
[38| to appear as a diamagnetic color current due to a re- 
sponse to an effective Yang-Mills magnetic field produced 
by SO interactions, which provides an explicit realization 
of a non-Abelian Landau diamagnetism. The equilibrium 
spin current in a 2D electron gas with a slightly modu- 
lated Rashba parameter was shown (39| to transfer spin 
from areas where spin is produced to areas where spin is 
absorbed. It was recently shown [4oT | that an equilibrium 
spin current in a 2D system with Rashba SO interactions 
results in a mechanical torque on a substrate near edge 
of the medium, which provides an experimental tool to 
detect the equilibrium spin current. Therefore, it can be 
concluded that a relation of the equilibrium spin current 
to spin transport should not be ruled out. 

The central result of the paper is an appearance of out- 
of-plane equilibrium spin current in the quasi-2D elec- 
tron gas under in-plane magnetic field in the presence 
of the SO interactions. In the absence of the magnetic 
field, the average values of the spin currents J Sai and 
3 s " are shown to coincide with the well-known results 
[HIH, El SI obtained for a strictly 2D electron gas, 
revealing a cubic dependence on the SO coupling con- 
stants. We show that the magnetic field strongly changes 
the in-plane spin-current components, contributing new 
terms to them. The new contributions turn to be pro- 
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portional, in addition to the magnetic field, either to the 
gate electric field or to Zeeman splitting. The new terms 
depend linearly on the SO coupling constants in the lim- 
iting case if one of the SO coupling constant is vanish- 
ingly small. Therefore, these contributions may prevail 
over the cubic dependent on the SO coupling constants 
terms. The out-of-plane component J Sy — {0, 0, J z '} 
vanishes completely with the magnetic field, and depends 
quadratically or linearly on the SO coupling constants. 

The paper is organized as follows. In Section II of 
this work we describe explicitly an analytical solution 
of quantum mechanical problem of one particle, moving 
in a quasi-2D system with finite thickness under in-plane 
magnetic field and in the presence of Rashba and Dressel- 
haus SO interactions by imposing a parabolic confining 
potential in the transverse direction. We take into ac- 
count in this work a gate potential too, which produces 
the SIA and Rashba SO interaction. In Section III we 
calculate the spin current in equilibrium. Conclusions 
are given in Section IV. In Appendix we present some 
routine calculations of the persistent spin current. 

II. ENERGY SPECTRUM IN THE PRESENCE 
OF AN IN-PLANE MAGNETIC FIELD 

We consider a quasi-2D gas of electrons, moving un- 
der an external in-plane magnetic field in the presence of 
both Rashba and Dresselhaus SO interactions and a gate 
potential. Single particle Hamiltonian of the system in 
the effective mass approximation can be written as 

~ !P ttl* uj 2 z 2 ~ I 

H= 2^ + ~^t~ ~ eEgZ + Hso + 2 mB<TB (1) 

where P = p— - A is an electron momentum in the pres- 
ence of a vector-potential A, m* and e are the electronic 
effective mass and charge, respectively; E g is a strength 
of the gate electric field. The second term in Eq. (JXJ) is 
the confining potential in z-direction, approximated as a 
parabola with a frequency ujq, which is a characteristic 
parameter of the electron gas thickness. This potential 
does not produce a structural inversion asymmetry and, 
consequently, SO interaction. Since Rashba SO interac- 
tion in the conduction band of a semiconductor is deter- 
mined by the electric field in the valence band rather than 
by that in the conduction band 43], the parabolic con- 
finement approximation neglects a small interface contri- 
bution to Rashba SO coupling constant. The last term in 
Eq. (|T|) is Zeeman splitting energy in the external mag- 
netic field B with u) z h = gfisB/2, where [Ib = 2^ i s 
the Bohr magneton of a free electron with mass mo , g is 
the effective Lande factor, and a — {cr x ,a y ,a z } are the 
Pauli spin matrices. 

Spin-orbital Hamiltonian H so in Eq. ([1} contains 
Rashba 0,3] term, Hr, due to a macroscopic SIA and 



Dresselhaus term [B|, [l5| , Ho , due to a BIA in the crys- 
talline structure. Dresselhaus SO interaction in bulk 
semiconductors with a zinc-blende crystal symmetry is 
proportional to the third order of the electron momen- 
tum P 

= \ E °MPhi - P ?+2)> (i = x,y,z;i + 3^ i), 

i 

(2) 

where 77 is a characteristic bulk coefficient of the SO split- 
ting. Since the average wave vector in the direction of the 
quantum confinement z is large, the terms involving p 2 
will dominate in Dresselhaus SO coupling for a quasi- 
2D electron gas with finite thickness. The expression for 
SO interaction Hamiltonian in MOSFETs and quantum 
wells of the width d grown along [001] crystallographic 
axis reads 

H so = H R + H D = -(<7 X Py ~ (TyP x ) + -{a X P X ~ O-yPy), 

(3) 

where a and j3 — —r)(p 2 ) — — j](ir/d) 2 are the sample de- 
pendent parameters of Rashba- and Dresselhaus-SO in- 
teractions, correspondingly. Rashba coefficient a is pro- 
portional to the gate electric field. Spin-orbital interac- 
tion can be interpreted as an interaction of a spin with 
randomly oriented, in accordance with the electron wave 
vector, effective magnetic field, which lies in the plane of 
the electron gas: 

H R = \o-n* n with Of // = |(Pxz), (4) 

h D 

H D = -a-U e f f , with 

= %{p^-{p^)),p y m)-p^M- (5) 

Although the effective magnetic field Q^ff, corre- 
sponding to Rashba term is perpendicular to the 2D wave 
vector of an electron, x-component of is in the same 
direction as p x while its y-component is directed in the 
opposite to p y direction. 

An external magnetic field is chosen to be in the plane 
of the 2D electron gas, along x axis B = {B, 0, 0} under 
the gauge A = {0, —Bz, 0}. The magnetic field tends to 
polarize the electron spin in the x-direction due to the 
Zeeman effect, and creates an angular momentum due to 
the orbital motion. 

In order to solve Schrodinger equation ih 9 ^^^'*'^ = 
H^(x, y, z; t) with a spinor ^(x, y, z;t) = (^) one ex- 
presses the electron wave functions with spin-up 
and spin-down orientations as y, z; t) — 

e ikxX+ik y y^ ^ z ^^ wn i cn yields 
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ih 



ift 



dt 



h 2 d 2 

2m* dz 2 
h 2 d 2 
' 2m* dz 2 



+ 



+ 



h 2 k 2 1 

+ (kyhuiB - eE g )z + — - >ifo+ [{a + ifi)(k y + uj B m*z/h) + ((3 + ia)k x + huj z ]ipi;(6) 
■ + (kyftjjJB - c-E g )z + ^— ^ + [(a - ifi)(ky + ujBm*z/h) + ((3 - ia)k x + hw z ]^,(7) 



r 



where ljb = eB/m*c is the cyclotron frequency, lo 



■ lUq is the effective frequency, and k = ^/k 2 



is the modulus of a 2D-wave vector. In the station- 
ary case the wave function is chosen as ^(x, y, z; t) — 
exp(—iEt/h)ty(x, y, z), where E is the total energy of an 
electron. In the absence of SO interactions and Zeeman 
term, Eqs. ([6]) and (J7J are decoupled and reduced to the 
oscillator equation with real wave function 



(z-z ) 2 -iH^ )(z) = 0, (8) 



h 2 d 2 

2m* dz 2 



where zq is the z-coordinate of a magnetic orbit, zq 



eE a —k v hu£ 



and E = E 



h 2 k 2 



(eE g -hkyU B ) 



is the en- 



ergy spectrum of the quantized orbital, E n = tuv(n+l/2), 
corresponding to the nth state 



Z~ Zq 

as 



(9) 



with H n (z) and = ^h/m*ui being the Hermite poly- 
nomial and the Bohr radius, correspondingly. 

SO interactions in Eqs. (|6|) and ([7]) mix the transverse- 
quantized levels, yielding the complex wave functions 



V>f-(z) and ipl( z )i they furthermore satisfy the condition 
V-"t4 = e y\ f ^ coordinate-dependent term, cx LOsm* z, 
in the off-diagonal part of Eqs. ^ and (J7J is origi- 
nated from the orbital magnetic field effect, which links 
nth orbital of a spin-up electron with (n ± l)th or- 
bital of a spin-down electron and vice versa. Eqs. ([6]) 
and (J7J are easily solved in the absence of this spatial- 
dependent term. Indeed, let us replace z in the 'mix- 
ing' terms of Eqs. © and ([7]) by the coordinate of 
the magnetic orbital center zq, and seek the solution as 
^ ( n\x, yi z) = e? k * x+ik y y i}£\z)(£). We get the follow- 
ing system of equations for A n and B n 



(E/Hlu -n - 1/2) An - c B n = 
c*A n + {E/hu -n- l/2)B n = 0, 



where the dimensionless coefficient cq is given 



(10) 



co 



(a + i/3) ( fcj,-| + ' 



huo 2 



+ (ia + f3)k x + u> z h 
(11)' 

The energy spectrum is immediately obtained from 
Eq.CGI 



E^(k x , k y 



2m* 



2m*u! 2 
{kyUi 2 + eE g w B /h) 



+ coh(n + 1/2) ±{ (a 2 + /3 2 ) 



h 2 k 2 (kyhujB -eE g f 



k 2 , {Jv4 + eE g uj B /K) 2 



(kyUj 2 + eE g w B /h) 



Pk x 



1/2 



(12) 



r 



The coefficients A n and B n in the spinor are completely linear combinations of tpn^ (z) 
defined from the normalization condition |v4„| 2 + |-B n | 2 = 
1 and Eq. (TTUl) 



o-*or 



A n = ^= and B„=±4^. 
V2 V2co 



(13) ^i(z) = e 



n=0 VObVt^ 



E 



-.11 „ 1^— ^1; (14) 



n=0 sJa B ^2 n i 



as 

Z - Zq 

as 



> (15) 



where the coefficients a n and b n satisfy the normalization 
condition J2^ =0 (\a n \ 2 + \b n \ 2 ) = 1- From the condition 
General solutions of Eqs. (0 and (J7J are sought as that Eqs. (JB]) and ([7]) are complex conjugate each other, 
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one gets a n = e %e b* n and b n — e lt 'a* with 9 being a real 
phase shift. Therefore, J2 n \ a n\ 2 = S« \ b n\ 2 = §• It is 
easy to estimate the average value of the spin operator 
components S — {S x , S y , S z } — h/2{a x ,a y ,a z } over the 
stationary states given by Eqs. ([T4")) and (fT5|) 

(S x )k = | + b* n a n ), 



n 

(s,} k = §£(K| 2 -|U 2 ) = o, 



(S+) k = ^<&„, and (5~) k = hj^ (16) 

n n 

So a spin precesses around the normal to the plane, and 
(S z )k averages out to zero, whereas in-plane components 
of the spin take finite values. The average position of an 
electron in the confining potential (z) can be calculated 
by the same way, 

(*)k = -| 1 {( a n+l fl « + a n a n+l) + 

(b* n+1 b n + b* n b n+1 )} + z , (17) 

which means that an overlap between the neighbor- 
ing transverse-quantized levels shifts the center of the 
magnetic orbit of both spin-up and spin-down electrons 
equally, in addition to the magnetic- and gate electric 
fields shift zq, in z-direction. 

Equations for the coefficients a n and b n with n = 
0, 1,2,3,... can be obtained by putting Eqs. (|T4"f and 
(1T51) into Eqs. © and © 



1 

2 J 



- c 6„ - V2n ci 6 n _i - 
V2(n+1) ci 6„+i = 0, 

- CpOn - \/2n <2 n _l - 



V2(n+1) cj a„+i =0, (19) 



where Co is defined by Eq. (fTTj). and the coefficient ci is 
given as 



ci — (a + iff) 



2hw 




(20) 



Note that a n = b n = for n < in Eqs. (JlSj) and (|T9|1. 

It is easy to see that the approximate equations (fT0|) 
can be obtained from Eqs.(fT8]) and (|T9|) by neglecting all 
terms ~ c\. 

In the absence of the external magnetic field, B = 0, 
the expressions for Co and ci, given by Eqs. flTlT) and (I2U1) . 



are simplified 

co = 7~^[* a ( fc:c ~ + + and Cl = °' 

(21) 

and, as a result, a mixing between the transverse- 
quantized levels is left off (see, Eqs. (fT8]) and (HH)). A 
simple exact expression for the energy spectrum in the 
absence of the magnetic field is obtained 



2 1.2 



2m* 2m*ujQ 



hco (n + - 
±\J{a 2 + P 2 )k 2 + 4afJk x k 



(22) 



which is a particular form of Eq. (fT2"|) written at B = 0, 
since Eq. (|12[) is exact in this limit. 

By expressing b n in Eq. (fl9l) through a n ,a n ±i and 
substituting into Eq. (JTS]) we get an equation for the 
vector a = {ao, ai, a,2, ■ ■ ■ }• An equation for the vector 
b = {6qi fri) &2j • • • } is obtained by the same way; finally 
we get: 



Na = 
Mb = 



(23) 
(24) 



N and M are square penthadiagonal matrices of infinite 
order with non-zero entries N^j 7^ {Mij 7^ 0) only if 
\i — j\ < 2, and N = (M)*. Apart from the non-zero 
main diagonal N n ^ n , the matrix N contains the first two 
diagonals, N n>n ±i and N n ^ n ±2, above and below it, which 
are given as 



(18) " £ 




The energy spectrum has to be found from the secular 
equation, by equating the determinant of the matrix N 
to zero. The infinite penthadiagonal matrix is truncated 
down to the first n rows and n columns, the roots of 



which can be found by numeric methods, |28l . [2 ( 
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III. SPIN CURRENT IN EQUILIBRIUM 



charge current 



The equilibrium spin current in the previous activi- 
ties 37H42J has been studied for a pure 2D electron gas 
in the absence of an external magnetic field. This sec- 
tion is addressed to study the spin-current in a quasi-2D 
electron gas with finite thickness in the presence of an 
in-plane magnetic field and the gate potential. In order 
to write the continuity equation for the charge density 
p = e\^(z)\ 2 = iV'tl 2 + l^jj 2 an d for the 7-component 
of the spin density S 1 = (h/2)(^' a^), the Schrodinger 
equations and ([7]) are multiplied to their complex- 
conjugate components ^| or ip*, which yields 



(30) 



where J and J ^ are the charge- and the spin-current, 
correspondingly. A violation of the spin conservation in 
the system results in an additional source term (torque) 
G 7 [39j in the spin-balance equation. The components of 
the charge- and spin currents read 

eh 



J, 



2m* 



-* T (cr X lo)i* 



n m*c 



(31) 



and 



.H 1 = -I- [* t cr 7 V ? * - Vj^Vt*] 

3 Atyl* 1 1 J J 7 J 



Ct 4. 

-{*> 7 (o- x z )j + (<t x z ) j( j 7 ]*} + 

f {^[^7 + - ^T c M&<r-r*)- ( 32 ) 



where a — {cr x , — a y , 0}, z$ is a unit vector, normal to 
the electron gas plane, and Aj is jth component of the 
vector-potential. It is easy to see that the structure of 
the contributions coming from Rashba and Dresselhaus 
terms is similar to the Hamiltonian form given by Eqs. 
(|4]) and (0. The magnetic field has a contribution to 
the charge current as well as to the spin current, which 
consists with contribution to the charge current 44[ ob- 
tained by means of Hamilton method in the absence of 
the SO interactions. 

The equilibrium spin- and charge currents at T=0 are 
found by averaging Eqs. (|3"Tj) and (|32|) over the stationary 
states, given by the wave functions (fT^|) and (jT5j) . and by 
integrating over all occupied states, which yields for the 



n=0 J y ' 



d 2 k ea eB 



(2tt) 2 * h 
d 2 k (Hk x B + ia „ B-ia . , 

+ 1T a " bn + —n~ Kan > ' (33) 



e E 

n=0 

W = e (^E/ ^ k y+n {ax) -n {a y ) + u}B{z) 



71=0 ' 



E 



<i 2 /c ( hk y a + iB a — iB 



— r H r — <A> + — - — &*a„ + 



(2tt) 2 [to* 7i 



+ ws-zo + Wfifl5 V2(n+ 1) (a^ + i«n + &£+i&n) f , (34) 



(J,) = 



(35) 



Note that the equations (pJ3")) - (f3l)|) for the charge cur- 
rent components can be obtained according to (Ji) = 
e(vi) as well by using the Heisenberg equation of motion 

v< = % = 

The average values of the spin current components read 

as 



d 2 k , ./ hk x a a 8 



(2tt) 2 [ 2m 



a /3 



(36) 



W) = 



E 





n=0 


a 




+ 2 






" d 2 fc f 


-y 


(2^) 2 \ 






+ 2 


e Ji/ ' 2 6 


and 





rf 2 fc 

(2^ 



, ul eE g uj B 



771 771 U)* 



Si 



7/ fcy (jj^ j 

in * / , )2 



as 



\o// 



(37) 



^) = WE/(0^^TI)>< 



(27T) 



where j = x,y and (Jf*) = (jf y ) = 0, also (Jf*) = 
(Jy ) = 0; eij is a 2D antisymmetric tensor with com- 
ponents e xy — —e yx = 1. The expressions for (cr x )k,« m ., 
(<7y)k,n m , (cz)k,n m and their evident momentum depen- 
dences are calculated in Appendix. 

It is evident that the normal to the electron gas com- 
ponent of (J s -) in Eq. (155} arises exclusively due to the 
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transverse confinement and the in-plane magnetic field. 
Nevertheless this term does not accumulate a spin, an 
electron tunneling from nth to (n + l)th level is accompa- 
nied, according to Eq. (155)) , by reverse flow from (n+l)th 
to nth level. 

The spin-continuity equation p0[) contains the source- 
term G 7 due to a violation of the spin conservation, the 
components of which are given by the following expres- 
sions 

Gj = -i^{&[a x (z Q x V)k* - [(V x z ) x a]j&*} - 

B e 
- «|{* f (V x a)j^ - (V x <r)i* + B-^(A x cr)* + 

1 e 

+ -gfi B ^(B x cr)* + a-& ((z x A) x a) (39) 

2 c 



where A„ for n = 0, 1 indicates the sub-band index with 
Ao = — and Ai = +, and the sign ± shows the spin- 
branch index. In-plane momentum dependence of the 
energy spectrum, Eq. (|42p . is determined by the terms 
I c 1 2 and (cgci + cqc\) 2 



CO 



4a/3k x [ fcj,-| 



2uj z h 



k ^ 



(« 2 +/3 2 ) 
H 2 

hio 2 

eEgUj B 
hw 2 



k ^ 



eEgUJB 
hw 2 



(43) 



where V = {V a ,— Vj/,0} and cr = {cr x , — Oy, 0}. The 
averaging of x- and y-components of the torque over the 
quantum-mechanic states gives 



(G x ) = 0, and (G y ) = 0, 



(40) 



which is in consistence with a result in strictly 2D system 
39] in the absence of an external magnetic field. Nev- 
ertheless z-component of the torque is not averaged to 



z a uj B m*)(a* n b n + 



L * ^ ; , o, 1 r, a 

+o n a n ) - i{ak y + pk x + -g^ B B + — z uj B m ){a n b n 



-b* n a n ) - i^-a B uj B m* \/2(n + 1) (a* +1 6„ - b* n+1 a n ) + 



-a B u B m 



iV2(n + l) K +1 b n + b* n+1 a n )}. (41) 



CqCi + CqC-l 



2a/3k 




So, |c | 2 ~ 0(a 2 ,B 2 ), whereas (cgcx + c c*) 2 — 
0(a 4 , (3 4 ,a 2 B 2 ) in the absence of Zeeman splitting. 
Therefore, expansion of Eq. (|4"2"j) over small SO coupling 
constants up to quadratic in a, 8 terms yields 



4 0) 



h 2 k 2 (eE q - hw B k v ) 2 1 , 

2^- *W +2 to± N to ' (45) 



The integration at T = is taken over each momen- 
tum component i = x,y and for both spin branches in 
each transverse-quantized subband up to the Fermi level, 
— iCJi ± — — Kn m ±- In order to calculate the average 
values of the charge- and spin current components, given 
by Eqs. (l35 ]) - (p0| . evident expressions of a n and b n are re- 
quired. The Fermi level is assumed to be set between n m 
and n m + 1 subbands, so that all levels up to {K m , n m } 
are occupied with a n ^ 0, b n ^ for n < n m and 
in = b n — for n > n m . A simplest case, which takes 
into account the inter-subband mixing due to an interfer- 
ence between the SO interactions and in-plane magnetic 
field, is n m = 1. In this case equations (fT5j) and (fH?|) are 
simplified to the form, given by Eq. (|79[) in Appendix. 
The analytical expressions for the energy spectrum in the 
first and second transverse-quantized subbands are given 
according to Eqs. ([82] ) -([85 ]) 



h 2 k 2 

2m* 



{eE g 



2m*u> 2 



nw + A„ — x 



l + 4|c | 2 + 8 



ci| 2 T ^\c \ 2 + 2(c* c 1 +c c* 1 ) 2 ,(42) 



h 2 k 2 (eE Q - huj B k v ) 2 3^ „ 
2^- a 2m- U > + 2^ T|C "'^ (46) 



The limit of integration over the occupied states can be 
found by fixing the Fermi energy Ep in Eq. (1421) and 
solving this equation for the momentum. It is neces- 
sary to note that an interference between the gate electric 
field and the orbital effect of the in-plane magnetic field 
shifts the Fermi surface along k y axis (see, Eqs. (|45|) . 
(|46l) and (|43[) ). Furthermore, Zeeman splitting makes 
the energy spectra asymmetric along both k x and k y 
axes. Therefore, the integrations over k x and k y have to 
be taken, generally speaking, over asymmetric intervals 
-K x m ,± <k x < K% m>± and -K,^ ± < k y < K v n ^ ± . 

The integration limit is calculated in Appendix by fix- 
ing the Fermi energy and transforming the momentum 
components in Eqs. (|45p and (l4l)l) into polar coordinates, 
k x = k cos ip, k y — ksimp. We express Eq. (|110[) for k^ ± 

as ± 



k^ ± 5k, where 
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L0 2 k 2 

Sh- 



oos 2 ip 

,*2 



. 2 

■ sin (p 



eE„ujB . 



2m* 



[E F - ujh(n + 1/2)] (cos 2 tp + ^| sin 2 tp) 



+ -^r\( a2 +P 2 ) ( cos2 <P + -T sin2 <A +4a^^sin< < 9Cos^l 
nr \ cj 4 / w 2 



1/2 



to 



ft 2 (cos 2 tp + ^ sin 2 (/J) 



oj 4 \ w 2 

(a 2 + /3 2 ) I cos 2 w H j sin 2 y> I + 4a/3— ^ sin y> cos 99. 

' w 4 / UJ Z 



(47) 
(48) 



for n = 0, 1. 

The average values of the Pauli spin-matrices 
(oj)k,n m i j — x i Vi z are calculated in Appendix. By using 
the expressions JIT), (J20J), flUSJ) and in Eqs. ([SI]), ® 
and ()98[) in Appendix, one gets an explicit momentum- 
dependent expressions for the averaged Pauli matrices 



( <J y)k±,r, 
(o"z)fe ± ,i 



= ±A„ 



ak + b 



Vck 2 + dk + e 

ak + b 
y/ck 2 + dk + e 

ak + b 
\Jck 2 + dk + e 



(49) 
(50) 
(51) 



where k is the momentum modulus {k x ,k y } = 
{k cos </>, k sin </>} in spherical-polar system and 



a = 


a — ^ sin 4> + p cos 0, 

W 2 


(52) 


b = 


eEgWB , . 
a_ fiw 2 W 


(53) 


a = 


a cos + /3 — ^ sin 0, 
w 2 


(54) 


b = 


p hJ 1 ' 


(55) 


a = 


W B /m* 2 2 

—*V—*{<r-P)aa*<P, 
con V wft 


(56) 


b = 


V ~ 

w V wfi 


(57) 



2 



/ 2eE„io B 



sin + /3 cos j + ^ sin0 + a cos0j ,(58) 

a ( a^r sin + /3 cos 6 ) + /3 ( /3^r sin 
\ w 2 / \ a; 2 

(59) 



V sw 2 

+acos0^ + 2w 2 /i ^a— |- sini/) + /3cos0^ 



(60) 



In-plane magnetic field induces an inter-sub-band cou- 
pling terms, (o»)£. ^ with i = which give a con- 
tribution to the y-components of the spin current (|37l) . 
These terms are calculated in Appendix. By neglect- 
ing the small terms [cq 2 ci + Cqc\ + (c* -I- ci)|c | 2 ] and 
[cjfci - clc\ + {c\ - ci)|c | 2 ], Eqs. and (011) in Ap- 
pendix are approximated as 



la ) off ft 


* 1 

C* + Ci 

2(e-l) 


A„(c* + ci 


auj B 
uin 


/ TO* 










. c* - Cl 


= -iA„(c* 


*2(e-l) 



(61) 



(62) 



So, the inter-subband coupling terms (c x )^„ and 

(&y)i/i depend only on the sub-band index and do not 
depend on the spin-branch index. 



In order to calculate the equilibrium spin-current com- 
ponents, Eqs. (f3"fj)) - (|3"5)) are integrated firstly over the 
momentum modulus k by taking into account Eqs. (|49[1 - 
(|62|) . Routine calculations yield under this condition the 
following results for the equilibrium spin-current J Sa= = 
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<4 3 



*9 



PW - « ) + 



1 feE g u B \ 2 (l 



| a 2_ j9 2| (o 2 +j9 2 )] + 



7aV x |a 2 -/? 2 | 



a 2 +/3 2 / l« 2 - /? 2 | 



{i( a 2 + 2/3 2 -|a 2 -/3 2 |) + ^[(a 2 -/3 2 ) 2 



eEgixJB f 



{a 2 - /3 2 ) 

*2 



|a 2 -/3 2 | a 2 +/? 2 - |a 2 -/3 2 |' 



1 - 



(a 2 - /3 2 ) 



,4 _ o4- 



1 / eEgWB_ V 
487ra / 3 2 \ ^ 2 ft J 



|a 2 -/3 2 | a 2 + /3 2 - |a 2 -/3 2 | 



-3(a 2 -/3 2 )|a 2 -/3 2 |-6a 2 /3 2 l^-AI 



5a; 2 ft 2 
127ra 



(a 2 - /3 2 ) 



3a 4 



(63) 



(64) 



well as for J s » = { JP , Jy y , 0} 



m* 2 ajo 



M./"' - " 3 > + 5^ (^r ) { 3 «* 3 + 5 " 2 >i° 2 " " 2 i - <« 3 + ^ 3 >< 3 " 3 + 23 « 



8/3 2 (a 2 +/3 2 ) 



2^ 



(a 2 - /? 2 ) 



2aT 

' (a 2 -/3 2 ) 



247ra 2 /3 2 w 2 \ 

2 



2 i fl2 _ |„2 _ «2| 



67ra (a 2 — /? 2 ) 



37r?i 4 (u; 

(4 o/„2 



(a 2 -/3 2 )-^(a 2 + /3 2 ) 







/? 2 

j3m*uj 2 B 
irh 2 ujLUo 
\a 2 -/3 2 \ 



Ep ujH 



1 / eEgU) B \ ' 
487ra 2 /3 V ^ 2 h J 



3a 4 



5a 2 - 3/3 2 



3|a 2 -/3 2 | + 2(4a 2 -3/3 2 ) |a2 ^ 



(a 2 - (3 2 ) J 6tt/3 



1 - 



\a 2 -P 2 \ _ a 2 + f3 2 -\a 2 - /3 2 | 
(a 2 - /3 2 ) ^ 



(65) 



(66) 



The non-zero component of S z spin-current, J Sz = {0, 0, </f z }, induced by the in-plane magnetic field, can be presented 



m*cj B (a 2 -/3 2 ) 
47ra/3w?i 2 



(a 2 +/3 2 -|a 2 -/3 2 |) (^) + 



OLL) 7 h 1 



a 2 -/3 2 



TTLUoH 



(67) 



The above expressions for the persistent spin current are simplified considerably as the gate voltage and Zeeman 
splitting vanish , E g = uj z = 0, 



jS y 



m* 2 ujQ 



1 2>Tth A UJ 



m* 2 ujo 



{ 3 

J s * ={0,0,0}. 



n(n2 2N 771 Up 

W 3^T 
3^7 



■}/ ,2 

(a — p ) H 2~(a +p ) 

w 



(a 2 -/? 2 ) 



3a; 



B {a 2 +0 2 ) 



Trh 2 u>uJo 
f3m*w% 
irh 2 LuuJo 



E F -~uh), 



(68) 



E F --uh), } ; (69) 



(70) 



Although the spin-current components turn to vanish at 
a = (3 in the absence of the magnetic field, the latter 
destroys this symmetry and yields a finite spin current 
as a — > (3. Indeed, the expressions (f63"l) - (|67[) yield the 
well-known results of Refs. [U [42[ for spin current in a 



2D electron gas in the absence of the magnetic field 

{*2 *2 *l 

J s - ={0, 0, 0}.(73) 
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It is easy to check that the diagonal components of the for arbitrary a/0. Nevertheless 
spin current vanish, (J^") — (Jy y ) = (Jf M ) = as B — > 



<J*> = I 1 



3w| 



3nh 4 L 
m* 2 ujoa 3 



am*ui B f 3 \ t 29a ( eE g ujB 
irh 2 ujuj a \ F 2 



24tt I Huj 2 



11a / eE„ujB 



(74) 
(75) 



3irh 4 uj 24tt 

under these conditions. On the other hand the diagonal spin-current components are nonzero for a — > and 8 ^ 




m* 2 w f3 3 t 11(3 (eE g uj B 
247 



37rS 4 o; 

m* 2 wo/3 3 
3tt/i 4 w 



1 



huj 2 
3w| 



2nujoh 



0; , 



8m*0J B 



Ep w/i 



/ eEgUj B 
3tt \ Hu>q 



; oV, 



(76) 

(77) 
(78) 



r 



New contributions to the pure 2D (~ a 3 ,/3 3 ) spin- 
current expressions in Eqs. (|63| -([66 |) are caused by 
electron-transfer mechanism between nearest-neighbor 
transverse-quantized sub-bands, induced by the in-plane 
magnetic field, and they vanish, consequently, with mag- 
netic field. These terms are proportional either to the 
gate electric field or to Zeeman splitting too. Indeed, the 
gate electric field changes, on the one hand, the electronic 
energy spectrum, and shifts, on the other hand, the cen- 
ter of an magnetic orbit zq along z-axis. The coefficient 
Co, which has a physical meaning of probability amplitude 
for an electron transition from one spin-polarized branch 
to other one in the same transverse-quantized subband 
n, according to Eqs. (fT8")) and (|19p . parametrically de- 
pends on Eg and lo z . The new terms in the spin-current 
expressions depend linearly on the SO coupling constants 
in the simplest case if one of the SO coupling constant 
is zero. Therefore, these terms may dominate over the 
pure 2D terms for some values of a, /3, E g and magnetic 
field. 



IV. CONCLUSIONS 

In this paper the equilibrium spin current is calculated 
for a quasi-2D electron gas with finite thickness under 
in-plane magnetic field in the presence of Rashba- and 
Dresselhaus spin-orbit interactions. Note that the prob- 



37J-|42| 



lem has been solved for 2D electron gas in Refs. 
the absence of the magnetic field. Our calculations show 
that the in-plane magnetic field generates out-of-plane 
spin current, which appears exclusively due to sub-band 



mixing by means of in-plane magnetic field. Although 
the equilibrium spin current vanishes at a = B in the ab- 
sence of the magnetic field (see, Refs. [4l|, H§|), in- plane 
magnetic field destroys this symmetry, yielding non-zero 
persistent spin-current at a — B. The magnetic field 
strongly changes 3 S:c and 3 Sy spin-current components, 
and contributes new terms to them. 
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APPENDIX 

Equations (fl8|) and (fT9|) are reduced for n = 1 into the 
following form 

(e - 3/2) oi - Co h - V2ci b = 
(e — 3/2) b\ — Cq a\ — \2c\ ag = 
(e - 1/2) a - cq b - V2 Cl h = 
(e - 1/2) b -c*a - V2cl a x = 0, (79) 

where e = E/(h uj), and the higher order terms 02 and 
&2 are neglected, since the levels n = and n = 1 only 
arc filled. Excluding a\ and b\ from these equations one 
gets 
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Eqs. ([79")) yield a system of equations for a\ and 61 too, which differs from Eqs. ([50)1 and (|5Tj) by interchanging 
the coefficients in the front of 00 in Eq. (|80|) and of 60 in Eq. (|8Tj) . Both system of equations results in the following 
expression for the energy spectrum 

[(e - 3/2)(e- 1/2) - 2| Cl | 2 ] 2 - |c | 2 [(e - 3/2) 2 + (e- 1/2) 2 ] - 2(cS 2 c 2 + c 2 cf) + |c | 4 = 0. (82) 

I 



By using the relation (e - 3/2) 2 + (e - 1/2) 2 = 2(e - 
3/2) (e- 1/2) + 1 in the second term of Eq. (JS2| one gets 
a quadratic equation for z = (g — 3/2) (g — 1/2) 

z 2 - 2(|c | 2 + 2| Cl | 2 ) 2 - |c | 2 + |c | 4 + 4| Cl | 4 - 
-2(cS 2 c 2 + clef) = (83) 

with the following solutions 

z = |c | 2 + 2| Cl | 2 ±^J\c \ 2 + 2(c* Cl + c c*) 2 . (84) 

Solution of Eq. (|84[l for the dimensionless energy g is 
expressed as 



where A n = ± indicates the first (Ao = — ) and second 
(Ai = +) energy subbands, and =F assigns two spin- 
polarized branches in each energy subband. Note that 
the coefficients Co and c\ are given by Eqs. (fTTj) and ([20]), 
furthermore cq only depends on the in-plane momentum 
components {k x ,k y }. 



A rl / I The coefficients a n and b n are complex parameters. 

g ± =1+yV 1 + 4 ' C °' 2+8 ' Ci ' 2=f4 V M 2 + 2(cSci +c c|) 2 , Solutions of Eqs. (JHOJ and dHB) with the normalization 

(85) condition yield for the modulus |ao| = \bo\ 



\ a o\ 



(iz IKiiMiz l)( g - 5) - i c °' 2 ~ 2 ' ci ' 2 ] 2 ~ W- IKiz 5) - 2 ' Ci ' 2 ] 

4(?- !)[(?-!)(?- |) - Icol 2 - 2| Cl | 2 ] [(?- §)(?- |) - (e- i)(|c | 2 + 2| Cl p) + |co| 



(86) 



In the absence of the magnetic field B 4 or ci ^ n = 1 can be defined from the relations (|80|) and (f8Tj) 
inter-subband coupling disappears, and Eq. (f84| is 
reduced to the form 



[(?-3/2) 2 -|c | 2 ] [?-l/2] 2 -|c | 2 ] =0 



(87) 



The expression for |ao| 2 is simplified as |ao| 2 = 1/2 under 
this condition. 

The symmetry relations a n = e 6* and b n = e lS a* for 
the complex coefficients a n — |a ?l le*^" and b n — \b n \e 1 ^ 
imply that a total phase 8 = + <fi b n is undefined param- 
eter, whereas the relative phase 0° — 0„ for n = and 



^ = ±exp{i(^-<)}, n = 0, 1; 



arg 



£o 1 

rf V 2^ ci V 2 ; 



(g-3/2)(g-l/2) + cf _ 



ci 2 ; c? 1 2^ 



(g-3/2)(g-l/2) + cf _ 



(89) 



(90) 
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The wave function can be expressed as 

y, z) = exp{i(k x x + k y y) — (z — z ) 2 /2a 2 B } 
H n ((z - z )/a B ) ( I 



E 



, (91) 



where the signs ± correspond to two different spin- 
polarized branches. 



In order to estimate the equilibrium charge- and spin- 
currents we have to calculate the mean values of the Pauli 
matrices (c7i)k,n m m the eigenstates given by Eqs. (TT4")) 
and ([T5j) up to the Fermi level n = n m , and integrate the 
results over {k x , k y } components of the in-plane momen- 
tum vector, —K\ i ± < ki < ±. Routine calcula- 
tions yield the following results for ((Jx}u.,n m 



and for (a y }k^ 



(o-x)k,n m = E ( a n b n)Vx\ 
n=0 



— a^bo + b^ao + a\bi + b\a\ = 



(4 + cp) [(e |) 2 + (g - \f - 2|c | 2 ] + 4(cgc 2 + c cf) 
4(e- l)[(i- |)(e- |) - |coP - 2| Cl p] 



(92) 



n=0 



= -i(ag6 - b^ao + a\b\ - b\ai) 



(c* - c )[(g- |) 2 + (g- I) 2 - 2|c | 2 ] - 4(cgc 2 - c cf ) 
4(g-l)[(g-|)(g-i)-|c | 2 -2| Cl p] 



(93) 



By neglecting the small terms 4(cqC 2 +coc 1 2 ) and 4(cqC 2 — 
cqcI 2 ) in the numerators of Eqs. AMI) and (|93p . corre- 
spondingly, and by using the expressions (|83[) and (j85|) 
for the energy spectrum, one gets for (a x )k ± ,n 

(a x ) k± . n = ±A„4f^^ (94) 
2|c | 

and for (a y )k ± , n 

K) fe± ,„ = TKi^f^- (95) 



y-componcnts of the spin-current contain in addition 
an inter-subband coupling terms (c^)£ £ 



E VW^K + iK + i)^( a b n )=^lbo + Ka ) = 

(cj + ci) [(g - I) (g - 1) - 2| Cl | 2 ] + (cg 2 Cl + c 2 ct) 
2(?-l)[(?-|)(^^)-| C o| 2 -2| Cl | 2 ] ' 

(°v)£L = E v/^+TTK+i^+r)^ " = -i«bo b* iao ) = 

■ ( c i - gl) [(g ~ I) (f ~ |) ~ 2| C1 1 2 ] + (cS 2 d - C 2 C» 
1 2(g-l)[(g-f)(g-i)-|c | 2 -2| Cl | 2 ] • 
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The transverse component of the spin current (J Sz ) is proportional to (<T z )k iTSm , which is given by a simple expression 
K)k,n m = V / 2(n+T)(< +1 &; +1 )o- 2 I n ) = V2(a iao - 6*fe ) = 



(CpCl -c c|) _ cgci - c ct f . 

2(e- 1) [(e- i) - | Co p - 2| Cl p] A » | Co | ' W 



I 

The momentum dependent factor CqCi — cqc\ in (<r z )k,n m into polar coordinates fc x = fccos^, k y = fcsintyS and 
is given as write Eqs. (|45|) and (j46]) with fixed Fermi energy as 

c^a-cocl = -i-^ x [™L[( a 2 -p 2 )k x -fa z ] . (99) (« 2 - A)) 2 = A x k 2 + A 2 « + A 3 , (100) 



(coh) 2 V wfi 

i the limit o 

space we transform the momentum components {k Xl k y } mentum modulus, and 



In order to find the limit of integration in momentum where k — k -\ h ^ B cog2 ip+ (jiJ' UJ 2 ) sin 2 <p ^ s a shifted mo- 



A = 

A 1 = 
A 2 



[E F - tf (n + 1)] (cos 2 y + g sin 2 g) + ^ 

^(cos 2 ^+^sin 2 ^) 2 
(a cos <p + siniyj) 2 + (/?cos</? + sin</?) 2 



4 m 



W(cos 2 ip + £| sin 2 </j) 2 



(101) 
(102) 



8eE g m* 2 uiB / '(H 5 uj 2 ) cos<p 
(cos 2 </? + % sin 2 (/J) 3 

2 

8w z m* 2 (a sin y> + j3 cos (/?) 
ft 4 (cos 2 ^+^sin 2 <^) 2 



( uj 2 uj 2 1 

<^ 2a/?(cos 2 </? ^sin 2 yj) - (a 2 + /3 2 )— | sin t^cos ip \ 

y lo- 1 bj 1 J 



(103) 



A, = 



( ^ ) cos 2 (p(a 2 +/3 2 -4a/3sm^cos^) 8 eEg m* 2 u B Uz cos<p(acos^ - /3sin^) 



(cos 2 </5 + ^ sin 2 y)) 4 



/i 5 u; 2 



(cos 2 + ^| sin 2 i^) 4 



4m* 2 a^ 



?i 4 (cos 2 ^+ ^|sin 2 <^) 2 



(104) 



The parameter Ap does not depend on small SO coupling constants a and f3 and weakly depends on the gate 
electric field and in-plane magnetic field, whereas A\ ~ 0(a 2 ,f3 2 ), A-i ~ 0(E g ujBCt 2 , E g u!B[3 2 ), and A3 ~ 
0(E 2 u 2 3 a 2 , E 2 lj 2 3 2 ). By introducing an unknowing parameter t/ Eq. (jlOOl) can be rewritten as 



(k 2 - A> + y) 2 = (Ax + 2y)n 2 + A 2 n + A 3 - 2yA + y 2 = {Ax + 2y) 

The parameter y is found under the condition 

R( y ) = or 

y 

yielding 



A? 



2(Ax + 2y) 



+ R(y). (105) 



3 Ax - 4A 2 AsAx-AAl 
1 2/ - (A Ai - A 3 )y + 



0, 



y 



4A - Aj 
6 



Q ± F* 

2 V 4 27 



P 



(106) 



(107) 



2 ^ V 4 ^27 
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where 



P = -A A 1 + A 3 



(4A - At 
12 



(AAp-A.r {A A l -A,){AA -A 1 ) AA Z A X ~A\ 
Q - 108 + 6 + 8 • (108) 



Expression for y can be simplified to the form 



y = 2A ° - 2^- (109) 



Finally, the Fermi momentum k„± for each sub-band and spin- branch is found from Eq. (|105[) under the condition 
R(y) = 



eE g uj B sin ^ , I a , A i A 3 , / A i , ^3 , A 2 



n,± fuv 2 cos 2 (^+4sinV V 4 4A ° V 4 4A ° 4 ^° 



The equation (|110[) with Eqs. (|104l) yields an evident 
expression for ki[±. 
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